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ABSTRACT 

A phenomenological turbulence model for kinetic Alfven waves in a magnetized collisionless 
plasma, able to reproduce the non-universal power-law spectra observed at the sub-ion scales 
in the solar wind and the terrestrial magnetosphere, is presented. The process of tempera¬ 
ture homogenization along distorted magnetic field lines, induced by Landau damping, affects 
the turbulence transfer time and results in a steepening of the sub-ion power-law spectrum of 
critically-balanced turbulence, whose exponent is sensitive to the ratio between the Alfven wave 
period and the nonlinear timescale. Transition from large-scale weak turbulence to smaller scale 
strong turbulence is captured and non local interactions, relevant in the case of steep spectra, 
are accounted for. 


Subject headings: plasmas — turbulence — waves — magnetic fields — solar wind 

1. Introduction 


Spacecraft measurements bo th in the solar wind 
and the Earth magnetosphere (iBruno fc Carbone! 


2OI3I : Alexandrova et all2013 . 2008a|l show power- 

law energy spectra for magnetic turbulent fluc¬ 
tuations. At MHD scales, where kinetic ef¬ 
fects are subdominant, observations support an 
Alfvenic energy cascade where the magnetic fluc¬ 
tuations transverse to the ambient field display 
a spectrum close to the prediction based 

on a “critical balance” (IGol dreich fc Shridhar 
1995 : Nazarenko &: SchekochihinT 2011 1 between 


the characteristic times of the nonlinear trans¬ 
verse dynamics and of the Alfven wave propaga¬ 
tion along the magnetic field lines. At sub-ionic 
scales, a power law is also observed in a range ex¬ 
tending from th e proton {pA to the e lectron {pA 
Larmor radius (ISahraoui et al. 20091 2010 . 2011 


lAlexandrova et alT 20121 : Chen et al. 2013 1. but 

the exponent appears to be less universal, with 
a distribution peaked near —2.8 a nd covering the 
interval [—3.1,—2.5] (Fig. 5 of ISahraoui et al 


(|2013f) b For comparison, the magnetic spec¬ 
tral exponent i n the MHD range i s estimated as 
— 1.63 ± 0.14 bv ISmith et al. ( 2006ll . 

Gyrokinetic simulations at /3 = 1 display a sub¬ 
ion power-law spectrum wit h comp arable expo- 
nents (~2 -8 in Howes et HI ( 2011bll or —3.1 in 
iTold et n! ( 2015 11. Fully kinetic par ticle-in-cell 


(PIC ) codes with a reduced mass ratio (IWan et al 


l201,5l) . a s well as hybrid Eule rian Vlasov-Maxwell 
models ( Servidio et al. 20151) also predict steep 
spectra at small scales, associated with coherent 
structures and deformation of the particle distri¬ 
bution functions. 

At sub-ion scales, two types of waves play a dy¬ 
namical role: whistlers for which ions are approx¬ 
imately cold and kinetic compressibility is negli¬ 
gible, and (low-frequency) kinetic Alfven waves 
(KAWs) for which density fluctuations are signif¬ 
icant. Reduced fluid-like models for the nonlinear 
dynamics of such waves have been developed and 
numerically simulated, leading to a —8/3 sub¬ 


ion spectral expo n ent (ISchekochihin et a 
iBoldvrev fc Pere'3 2012: iBoldvrev et al. 


2009; 


2013: 


1 




















































































Mevrand fc Galti^ l2013[l . Nevertheless a com¬ 


prehensive understanding of turbulence at these 
scales is still missing. Carrying the critical-balance 
phenomenology to the KAW cascade leads to a 
' energy spectrum, significantly shallower 
than observed. Still assuming critical balance, 
Boldyrev et al.1 ( 2013l l proposed a steeper spec¬ 


trum resulting from coherent struct ures and in- 
termit tency corrections. Differently, iHowes et al 


( 2008ll suggested a balance between Landau damp¬ 
ing and energy transfer. The spectrum is 

then multiplied by an exponential factor originat¬ 
ing from the variation of the energy f lux along the 


cascade. This model was extended in Howes et al 


(l 2011 ar) . in an attempt to include nonlocal inter¬ 
actions, relevant for steep spectra. 

In this letter we revisit KAW phenomenology, 
including in the energy transfer time the effect of 
ion temperature homogenization along magnetic 
field lines induced by Landau damping. The lat¬ 
ter, whose linear rate i s assumed to persis t in the 
nonlinear regime (see ISchekochihin et Hi ( 201 ^ 
for a discussion), not only provides the exponential 
cut-off (at a scale depe nding on the Kolmo gorov 
constant) discussed in iHowes et al.l (l2006ll . but 
also leads to a steepening of the —7/3 power law. 
Furthermore, as it does not necessarily enforce 
strict critical balance, this model describes the 
transition between weak turbulence at large scales 
to strong turbulence at smaller scales. 


2. Model setting 

We consider a collisionless proton-electron 
plasma permeated by a strong ambient magnetic 
field (of amplitude Hq), with equal and isotropic 
mean temperatures Tg = Ti. Alfven waves are 
driven in the MHD range, at scales much larger 
than Pi. The transverse magnetic field fluctuations 
SB± are measured in velocity units by defining 
b = va{SB±/Bo), where va is the Alfven velocity. 
The amplitude of these fluctuations at a trans¬ 
verse wavenumber k± is bk {kj_Eky^'^, where 
Ek is the transverse magnetic spectrum. 

2.1. Involved time scales 

The magnetic held being stretched by electron 
velocity gradients and the dynamics dominantly 
transverse, we dehne, when the interactions are 
local, a stretching frequency (inverse of the non¬ 


linear time t/vl) ujnl k±Vek where the trans¬ 
verse electron velocity Vek at scale kj^ is given by 
Vek = abk- Here, a is a function of k±pi, equal 
to 1 in the MHD range and scaling like kj_pi in 
the far sub-ion range, with a smooth transition 
near the ion gyroscale. As in hydrod ynamic tur¬ 
bulence, we write lun t, (iKovasznay 

1948 : Pancheylll969t i. Neyertheless, when Ek is de¬ 


caying fast enough, the above expression does not 
necessarily ensure the expected monotonic growth 
oi ujnl- In this case nonlocal interactions cannot 
be neglected, and ujnl should rather be viewed 
as the stretching rate due to all the scales larger 
than kj^. It is then taken equal to the r.m.s. 
value ojnl = AJ E„ dp]^''‘^ (jElisson 1962; 


PanchevI Il971t iLesieun 1200811 , where A is a con¬ 
stant. The local approximation is recovered when 
the integral diverges at large k±, while the in¬ 
tegral formula can be replaced by the equation 


duj%i^/dk± = A?c 


KEk. 


Alfven waves are characterized by a frequency 
uJw = wfeyWA and a dissipation rate 7 = jk^^VA- 
Here ZU (equal to 1 in the MHD range) and 7 
are functions of kj_pi provided by the kinetic the¬ 
ory. In a linear description, parallel and perpen¬ 
dicular wavenumbers are defined relatively to the 
ambient field (taken in the z-direction), but dis¬ 
tortions of the magnetic field lines should be re¬ 
tained in the nonlinear regime. We write fey = 

kz + p']_Epdp]'^/'^, where k^ is a constant 

estimating the parallel wavenumber of the energy- 
containing Alfven waves. It should not be con¬ 
fused with the inverse correlation length along 
the z-direction. Differently, measures the 

correlation along the magnetic field lines (which 
is larger than along the z-axis). A procedure 
to estimat e fey in numerical sim ulations is de¬ 
scribed in Cho fc LazarianI ( 2004 1. while an ap¬ 
proach us ing the frequency a s a p roxy for fcy 


is used in 


ng tne irequency a s a p rc 
TenBarge fc How^ ( 20121) . 


This sug¬ 


gests to write uiw = ojkzVA + dojw with a tur¬ 
bulent frequency shift dtow = [fo^ ^^P^Epdp]^^^ 
or d{dujw)‘^/dk± = uJ^k’^Ek- We similarly write 
7 = ikzVA + ^7 with d{6^)‘^/dk± = ^'^k']_Ek. A 
multiplicative constant should also enter the defi¬ 
nitions of 5uJw and 7 , but is easily scaled out. 

Another time scale originates from the com¬ 
pressible character of KAWs. The latter are sub¬ 
ject to Landau damping, resulting in tempera- 
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ture homogenization along the magnetic field lines, 
on the correlation length in a time Tur 
{vthrk\\)~^■ Here the thermal velocity Vthr ap¬ 
pears as the r.m.s. streaming velocity of the r- 
particles. This time scale, which arises explic¬ 
itly in Landau fluid closures (Hammett et al.lll997 : 
Snyder et al. 1997 : Sulem fc Passod 2015 ). is very 
short for the electrons, while for the ions it is com¬ 
parable to that of the other relevant processes and 
can thus affect the dynamics. Due to magnetic 
field distortion, this process introduces additional 
nonlinear couplings characterized by the frequency 
loh = where p, is a numerical constant. 

We are thus led to write ujh = IJ-Vthikz + Swh 
with d{5uiHY / dki_ = ji^pk^Ek- 

Finally, we define the transfer time Ttr = 
tnl {tnl/tw + tnl/th), or in frequency terms 
uJtr = -b w//). When only one pro¬ 

cess competes with n onlinear stretching, the criti- 
cal balance condition (iGoldreich fc Shridhaijll99,^ 


Nazarenko fc Schekochihinl2011tLSchekochihin et al 

20091 ) ensures the equality of the two associated 
time scales, and thus of the transfer and stretching 
times. 


2.2. The nonlocal model 


the small-scale eddies cannot be sufficiently highly 
correlated wit h each other to contribute equally 
( ElissonI 1962), we here use the local approxi¬ 
mation kj_Ek of the Reynolds stress rather than 
the original Obukhov’s expression Epdp which 
leads to an unphysical behavior in the dissipa¬ 
tion range of hydrodynamic turbulence (jPanchev 
19711) . 


Normalizing frequencies by Hi, wavenumbers 
by p~^, energy spectra by v\pi, energy fluxes by 
v\Vti, and denoting by /3 the ion beta, we obtain 
the non-dimensional model equations (keeping the 
same notations) 


duj%L/dk± = A^p-^a^klEk ( 1 ) 

d{SuJw)‘^ / dk± = f3~^Lj‘^k']_Ek (2) 

d{S^f/dk^ = r^tklEk (3) 

d{5ujH)'^ / dk± = p^k]_Ek (4) 

de/dkji_ =-2[l3~^/‘^^k:z + {S^)\Ek (5) 


Ek — [(/3 + p)kz -\- doJw “b duJH\-r -2—■ 

( 6 ) 


Except possibly near fc_L = 1, a = lo. The nonlin¬ 
earity parameter y = obeys 


Retaining KAW linear Landau d amping leads 


to the phen omenological equation (I Howes et al 


2008, 2011 


dtEk +Tk = -2^Ek + Sk, where 
Sk is the driving term acting at large scales and 
7/c the transfer term related to the energy flux e 
by 7fe = dejdk±. Due to Landau damping, en¬ 
ergy is not transferred conservatively along the 
cascade, making e scale-dependent. For a steady 
state and outside the injection range, one has 
de/dk± = —2jEk. Note that the present setting 
differs from the asymptot ic regime considered by 
ISchekochihin et ^ ( 2009h for which, under the si¬ 
multaneous conditions k±pi ^ 1 and k±pe <C 1, 
KAWs are not subject to Landau damping, as they 
transfer part of their energy via parallel phase mix¬ 
ing to the ion (electron) entropy cascades only at 
ion (electron) gyroscales, where it is cascaded both 
in physical and velocity spaces to collisional scales 
via perpendicular phase mixing. 


Estimate of the energy flux relies on a basic tur¬ 
bulence description (overlooking intermittency). 
We write e = CoJtrkj_Ek, where C is a negative 
power of the Kolmogorov constant. Arguing that 


dk± 


a^k^Ek 




f (7) 


where / = ^^(wW) + is positive. 

At the (small) injection wavenumber fco, turbu¬ 
lence is characterized by A = kz/lk^Eko]^^'^ = 
{kz/ko){Bo/dB±o) and, when taking = 

anda;(^) = p-^/^kz + k^^^E^^), 
Xo = A/{\ + A). In the strong turbulence regime 
with fcz = 0, X = A in the full spectral range, 
while for ^ 0, y starts growing near k^ but 
cannot exceed A. 


3. A simplified local-interaction model 
3.1. The usual conservative cascades 

Assuming local interactions, we have 

El ~ +kY‘Ei'\ 

( 8 ) 

When neglecting dissipation (constant e), we re¬ 
cover the usual inertial energy spectra. For weak 
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A 

0.71 

1 

1.22 

1.30 

1.41 

2 

4.47 

exponent 

-3.18 

-2.81 

-2.69 

-2.66 

-2.63 

-2.53 

-2.45 


Table 1: Sub-ion exponent versus A in conditions of Fig. [T^ (fitting range 8 < k± < 40). 


turbulence <C k, 

MHD range, while in the sub-ion range (as uJ ^ 
Ek oc For strong turbulence {kz 

negligible), Ek oc in the MHD range, while 

in the sub-ion range Ek(xkj_ .A kj_ regime is 
also obtained in the sub-ion range when the effect 
of wave propagation is negligible, as observed in 
two-dimensional hybrid PIC simulat ions with an 


_L < kz), Ek oc 


in the 


out-o f-plane ambient magnetic field (|Franci et al 
2 OI 5 II. Furthermore, as easilv seen from ea. ®, 


an additional regime with Ek oc is possible at 
very large scales when, for small enough e, turbu¬ 
lence is not yet developed and watl almost con¬ 
stant. Such a spectral exponent is observed in the 
solar w ind at scales larger than the k , inertial 
range ( Matthaeus &: Goldstein 19861 : Nicol et al 


200^IWicks et al.ll2ninll . 


3.2. Effect of Landau damping 

For strong turbulence, when assuming local in¬ 
teractions and neglecting kz contributions. 


Cfc = eo exp 


-2C'"^A“ 


I ko 


T' c- 
7=2 (w- 




(9) 

where the notation ek for the energy flux stresses 
its wavenumber dependence. This quantity is to 
be substituted in Ek - / 3 i/ 3 A- 4 / 3 c'- 2 / 3 g 2 / 3 ^- 7 / 3 ^ 
taking a = oJ and 7 /aJ^ ~ 6{f3) OJSpeJpi 
when using eqs. (62) and (63) of I Howes et al 


( 2006l l with /3 = 1). Furthermore, in this range 
tJ ~ a{P)k_L with a(/3) = (1 -I- /3)“^. This 
leads to ek ^ eokj^^ exp{—2a{f3)C~'^A~'^S{f3)kj_), 
with C = 25(/3)C'“^^A“^/3^/^. This results in 
a steepening of the algebraic prefactor of the 
magnetic spectrum which is now proportional 
to ;j“C/ 3+2C/3) exp[— (4/3)a(/3)(5(/3)C'~ ^A~^fc 1 ]. 
Comp ared with the spectrum obtained in iHowes et al 
(l2008ll . the present model predicts that, in addi¬ 
tion to an exponential cut-off, Landau damping 
leads to a correction of the power-law exponent. 

In contra st with intermitte n cy co rrections dis¬ 
cussed in iBoldvrev fc Perea (|2012ll . this correc- 



Fig. 1.— (a): Normalized energy spectrum Ek 
(black lines, l.h.s. labels and outer tickmarks) and 
parallel wavenumber k\\ (blue lines, r.h.s labels and 
inner tickmarks), for fcz=0,/3 = e= l and A = 
2 (dashed lines), 1 (solid lines) and 0.5 (dotted 
lines), (see other parameters in text); (b): Same 
as (a) for A = 1 and eq = 10“^. 
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tion is not universal, which is expected when dis¬ 
sipation and nonlinear transfer ti mes display the 


same wavenumber dependence (iBratanov et al 


2013ll . This situation holds for the second term 
in the exponential arising in eq. that leads 
to the correction of the —7/3 spectral index. 
We note that it defines a dissipation length 

oc and thus a 
dissipation rate Wd = Vekd for Cfc, which has the 
same fcj_ dependency as uinl that, by a critical 
balance argument, identihes with totr beyond the 
transition range. A similar power law decay of e 
is encountered in drift-kinetic plasma turb ulence 
(see eq. (2.53) of Schekochihin et al. ( 2015ll '). 

Differently, for weak turbulence, we get 


Cfc = 


- A-^C-^p-^k 


f'k± — 

f -L ^ 


'ko 




( 10 ) 

1/2 

where the integral behaves like fcj/ . Equation 
OT predicts that Ck and thus Ek vanish at a finite 
k ±, indicating the breaking of the analysis near the 
corresponding scale, an effect possibly related to 
the difficulty for weak turbulence to exist in the 
presence of a significant Landau damping. 

4. Numerical integration of the full model 

When the spectrum is too steep, a numeri¬ 
cal integration of differential equations (HD-® 
is needed. The functions uj and 7 are then 
evaluated from the full linear kin etic theory 
by m eans of the WHAMP software ( RbnnmarM 
19821 1. Moreover, a = oJ. We prescribed 
conditions at k± = ko in the form = 

57(0) = ^ ^^3/2^V2 

eo = eCkoio^^^j^ Eo/[{p-^/‘^uJ + ^)k:, + Sww + Swh]- 

Here Eq is an arbitrary constant (taken equal to 1 
with no lack of generality). We chose kg = 10“^, 
C = 1.25 and fj, = 1.8. Except when otherwise 
specified, we also took e = 1. For clarity’s sake, 
when several energy spectra are plotted in the 
same panel, one of them (solid line) is normalized 
by its value at fcj, = 1, while the others (dashed 
and doted lines) are rescaled to make all the spec¬ 
tra equal at k± = fcp. In red are indicated the 
fitting ranges (dashed-dotted straight lines) and 
the corresponding spectral exponents whose last 
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Fig. 2.— Same as Fig. la with P = 0.01, A = 1 
(solid lines)and 0.45 (dashed lines) (a), and P = 
10, A = 3.16 (solid lines) and 1 (dashed lines) (b). 
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digit only is sensitive to moderate changes of the 
fitting domain. 

In Fig. [T^, we focus on the strong turbulence 
regime for /3 = 1, assuming kz = 0- In the MHD 
range, a energy spectrum establishes and 

fc|| scales like fcj_ in all the cases. Differently, at 
the sub-ion scales, the spectrum is steeper when A 
is smaller, displaying exponents —2.53 for A = 2, 
—2.81 for A = 1, and a fast decay for A = 0.5. 
In this range, fey increases slower, and this even 
more so when A is smaller. In the present set¬ 
ting, A « 0.71 (leading to a spectral exponent 
—3.18) and A ft! 4.47 (exponent —2.45) appear as 
the extreme values for existence of an extended 
power-law spectrum at the sub-ion scales. These 
limit exponents are consistent with the dispersion 
of solar-wind measurements. Exponents for in¬ 
termediate values of A are displayed in Table 1. 
Their statistical distribution (and the most prob¬ 
able value) are nevertheless beyond the scope of 
the model. Note that a —7/3 exponent (rarely re¬ 
ported in observations) is approached for large A, 
only when ji — 0. 


When keeping A = 1 but decreasing the en¬ 
ergy transfer rate eo by taking e = 10“^ (Fig. 
[Ib), the MHD and sub-ion spectral exponents 
are not affected, but a range becomes vis¬ 
ible at the largest scales. In this range, where 
fc|| remains small, turbulence is not developed, 
possibly as in the solar wind energy containing 
range. Further decrease of cq leads to a kj^ 


range extending down to 
ation observed in the mae 

the ion gyroscale, a situ- 
netosheath near the bow 

shock( 

Czavkowska et al. 

2 OOII: [Alexandrova et al. 

2008bl. 


Another issue is the influence of /3, keeping 
kz = 0. For f) = 0.01 (Fig. [2b), we considered 
A = 1 and A = 0.45, leading to sub-ion spec¬ 
tral exponents —2.37 and —4.21 respectively, while 
as expected the MHD range is not affected. For 
/3 = 10 (Fig. Od), we used A = 3.16 and 1, for 
which the sub-ion exponents are —2.86 and —3.66 
respectively. For the former value of A, a spectral 
bump is visible in the transition zone, consequence 
of the sharp decrease of the ion Lan dau damping 


at the se scales (see e.g. Fig. 3 of iHowes et al 

ll^). 


To justify the chosen values of A (which is 
equal to x in the above setting where kz = 0 


and prescribes its saturated value when kz ^ 0), 
it should be noted that, at a fixed wavenumber, 
X increases with the amplitude of the fluctua¬ 
tions, linearly in weak turbulence and at a slower 
rate when the amplitude gets larger, a behav¬ 
ior supported by FLR-Landau fluid simulations 
(in preparation). Furthermore, the Mach num¬ 
ber, which scales like (where bo measures 

the amplitude of the large-scale turbulent fluctua¬ 
tions) is usually observed to be mod erate in the so¬ 
lar wind ( Bavassano fc Brunoll995ll . i n spite of the 


broad range of reported values of /3 (iChen et al 


20141 1. This implies that the turbulence level, and 


thus A, should be decreased at smaller /3. 

Figure [3] addresses the influence of the param¬ 
eter A when A is fixed at a moderate value, here 
•\/2. Increasing A results in a change from strong 
to weak turbulence near the driving scale. For 
A = 0.5 (Fig. |3b), the function y rapidly satu¬ 
rates to a value slightly smaller than A, establish¬ 
ing critical balance and thus a strong turbulence 
regime. The spectral exponent —1.72 measured in 
the MHD range differs from —5/3, since x is still in 
the growing phase, but choosing ko smaller would 
ensure a kj^^^ critically-balanced MHD range. At 
small scales, the exponent is —2.62, comparable to 
the values displayed in Fig. [Tb. The case A = 5 
corresponds to an intermediate regime where xo 
is significantly smaller, resulting in a —1.92 large- 
scale spectrum, steeper than —5/3 but neverthe¬ 
less distinguishable from the —2 weak-turbulence 
value. The sub-ion exponent —2.71 is consistent 
with a strong turbulence regime, characterized by 
an almost constant x at these scales. In Fig. [Sb, 
A = 10 and 15, both lead to a weak-turbulence 
regime at large scales, with exponents —1.95 and 
— 1.97, very close to the theoretical value. Dif¬ 
ferences nevertheless hold at small scales. For 
A = 10, the sub-ion dynamics corresponds to 
strong turbulence, qualitatively similar to the case 
A = 5, with a spectral exponent —2.9 and an al¬ 
most constant Xj although somewhat smaller. Dif¬ 
ferently, for A = 15, the fluctuations are too weak 
for the energy transfer to efficiently compete with 
Landau damping, leading to an exponential decay 
of sub-ion spectrum, and a function x which starts 
to decrease by k±pi = 10. 
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5. Conclusion 



Fig. 3.— (a): Energy spectrum Ek (l.h.s. labels) 
and nonlinear parameter x (r.h.s. labels) for A = 
•\/2 with A = 0.5 (solid lines) and A = 5 (dashed 
lines); (b): Same as (a), with A = 10 (solid lines) 
and A = 15 (dashed lines). 


This model predicts a non-universal power- 
law spectrum for strong turbulence at the sub¬ 
ion scales (beyond the transition range) with an 
exponent which, in contrast with the —5/3 in¬ 
ertial MHD cascade, depends on the saturation 
level of the nonlinearity parameter y, covering a 
range of values consistent with solar wind and 
magnetosheath observations. Such a non uni¬ 
versality, associated with Landau damping, was 
also reported in three-dimens ional PIC simula - 
tions of whistler turbulence ( Gary et al.l 120121 
and FLR-Landau fluid simulations of KAW tur¬ 
bulence (in preparation). As the present ap¬ 
proach does not capture sub-electron scale dy¬ 
namics, the exponential cut-off might be re¬ 
placed by another regime, such as the steep power 
l aws reported from both numerical simu l ations 
( Camnoreale fc BurgersI 2011 : Gary et al. 2012 ) 
and spacecraft data (|Sahraoui et al.1120131) . 


The research leading to these results has re¬ 
ceived funding from the European Commission’s 
Seventh Framework Programme (FP7/2007-2013) 
under the grant agreement SHOCK (proiect num¬ 
ber 284515). 
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